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Abstract 

We discuss a few, apparently different (but actually, tightly re- 
lated) problems: 

1. The relation between QCD and valence quark model, 

2. The evaluation of the nonlocal condensate {q{x)q{0)), its relation 
to heavy-light qQ quark system and to constituent quark mass, 

3. The asymptotic behavior of the nonperturbative pion wave func- 
tion ip{k'j_,x) at X — > 0, Ijk"^ —>■ oo and 

4. The large order behavior of perturbative series. 

The analysis is based on such general methods as dispersion rela- 
tions, duality and PCAC. We use the steepest descent method (also 
known as semiclassical, or instanton calculus), introduced by Lipatov 
to calculate the n— th moment of the ip{k'^,x) with result {k'^) ~ n\. 
This information fixes of the asymptotic behavior of wf at large k'j_. 
This behavior is turned out to be Gaussian commonly used in the 
phenomenological analyses. The same method determines the asymp- 
totic behavior of the mixed local vacuum condensates {qG^j^q) ~ n! at 
large n as well as the nonlocal vacuum condensate {q{x)q{^)) which 
naturally arises in the description of the heavy-light qQ quark system. 

The relation between nonlocal condensate and constituent quark 
mass is also discussed. 



^ On leave of absence from Budkcr Institute of Nuclear Physics, 
Novosibirsk, 630090, Russia. 
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1. Introduction 

The problem of bound states in the relativistic quantum field theory with 
large coupling constant is an extremely difficult problem. To understand 
the structure of the bound state is a very ambitious goal which assumes 
the solution of a whole spectrum of tightly connected problems, such as 
confinement, chiral symmetry breaking phenomenon, and many others which 
are greatly important in the low energy region. 

A less ambitious purpose is the study of the hadron wave function (wf) 
with a minimal number of constituents ( we consider such an analysis as the 
possible way of understanding the valence quark model in terms of QCD.). 
As is known such a function gives the parametrically leading contributions 
to hard exclusive processes. The corresponding wave functions within QCD 
have been introduced to the theory in the late seventies and early eighties 

to describe the exclusive processes. We refer to the review papers [0,|0|7 
0] on this subject for the details. 

The main idea of the approach is the separation of the large and 

small distance physics. At small distances we can use the standard perturba- 
tive expansion due to the asymptotic freedom and smallness of the coupling 
constant. All nontrivial, large distance physics is hidden into the nonpertur- 
bative wave function (wf) and can not be found by perturbative technique, 
but rather it should be extracted from elsewhere. The most powerful analyt- 
ical nonperturbative method for such problems, I believe, is the QCD sum 
rules §,§. 

The first application of QCD sum rules to the analysis of nonperturbative 
wf was considered more than decade ago 0. The information extracted for 
the few lowest moments, unambiguously shows the asymmetric form of the 
distribution amplitudes. Since then this subject is a very controversial issue 
||^-[|T^ and we are not going to comment these quite opposite points in the 
present note. However we would like to make a comment that a finite number 
of moments is not sufficient to completely determine the wf. The behavior 
of the asymptotically distant terms plays a crucial role in the reconstruction 
of the wf. 

• The main goals of the paper are: 

The analysis of the asymptotically distant terms at n — oo in the lon- 
gitudinal {{2x — 1)") and transverse {k"^) directions. We are going to use 
very general methods in this analysis, like dispersion relations and duality. 
Thus, the obtained results should be considered as the direct consequences 
of QCD. 

Besides that, we use the same methods to analyse the so-called nonlocal 
condensate (NLC) {q{x)P exp(ig Jq A^dXf^)q{0)) which is the fundamental 
object of the QCD-theory. It naturally appears in the analysis of the heavy- 
light quark system. 

We make contact between valence quark model and light cone wf defined 
exclusively in QCD terms. Finally, we define and calculate a constituent mass 
of the quark in terms of NLC, i.e. in terms of the object which is defined 
itself in terms of QCD. 

Let me emphasize from the very beginning that the ideology and methods 
(unitarity, dispersion relations, duality) we use are motivated by QCD sum 
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rules. However, we do not use the QCD sum rules in the common sense: we 
do not fit them to extract any information about lowest resonance (as people 
usually do in this approach), we do not use any numerical approximation or 
implicit assumption about higher states. Instead, we concentrate on analy- 
sis of the appropriate correlation functions themselves to extract the most 
general information. 

The paper is organized as follows. In the next section we define the 
nonperturbative wf through its moments. 

Section 3 is devoted to the analysis of the large order behavior. We 
use the steepest descent method (also known as semiclassical, or instanton 
calculus) introduced by Lipatov, to calculate the n— th order of perturbative 
series. After that the dispersion relations unambiguously fix the asymptotic 
behavior of w/ at large k'j_. 

In section 4 we apply the same method for consideration of the different 
object - the NLC, which originally was introduced for analysing the heavy- 
light quark system. 

2. The Definitions. The constituent quark model and QCD. 

We define the pion axial wave function in the following gauge-invariant 
way: 

tUq.Mzq, z') = (0|d(;^)7,75e^'^-^^^'^"n(-z)|7r(g)) (1) 
= E ^(0|J(0)7;.75(^^. 5.)"«(0)|7r(g)), 

where Dy=Di, — and iD^ = id^ + gA'^^^ is the covariant derivative. 
From its definition is clear that the set of different vr meson matrix elements 
defines the nonperturbative wave function. 

The most important part (at asymptotically high q^) is the one related to 
the longitudinal distribution. In this case 2;^ ^ and the wf depends only 
on one zq- variable. The corresponding Fourier transformed wave function 
will be denoted as 0(^) and its n—ih moment is given by the following local 
matrix element: 

(0M7.75(^ D, z,ru\Ti{q)) = thquizq^^C) = ifnquizqT d^Cm (2) 

— q^ ^ 00, ~ 1 ^ = xi — X2, xi + X2 = 1, z'^ = 0. 

Therefore, if we knew all matrix elements (|^) (which are perfectly well- 
defined) we could restore the whole wf. The QCD sum rules approach 
allows one to find the magnitudes only the few first moments 0. As is 
known, this information is not enough to reconstruct the wf; the parametric 
behavior at ^ — ±1 is the crucial issue in this reconstruction. To extract the 
corresponding information, let us consider the following correlation function: 

t J dxe''^^0\TJl{x), Jo(0)|0) = (zqy^+'lniq^), jl| = d^,z,j,{i D, z^)"^i(3) 
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and calculate its asymptotic behavior at large q^. The result can be presented 
in the form of the dispersion integral, whose spectral density is determined 
by the pure perturbative one- loop diagram (Fig.l): 

1 roo Iml^^rt^ /m/„(sr'-* = — ^- -. (4) 

nJo s-q^ ' ^' 47r(n + l)(n + 3) ^' 

Our main assumption corresponds to the very likely case that the vr meson 
gives a nonzero contribution to the dispersion integral for arbitrary n and, 
in particular, for n — > oo. Formally, it can be written in the following way 



/ dslmlisf^^"' = - / dslml{s)l, (5) 
vr JO 71 Jo 



We assume that 5" in this relation is some finite, nonzero number (we can say 
that " the pion fills" a finite duality interval). Our assumption means that 
there are no special cancelations and vr meson contribution to the dispersion 
integral is not zero, i.e. S^dl) 7^ 0, where we specified the notation for the 
longitudinal distribution. In this case at ^ 00 our assumption (|) leads 
to the following relation: 

/^{{"X"--)-^ (6) 
It unambiguously implies the following behavior at the end-point region 0: 

r) = d^cm ~ i/n', 0(e ^ ±i) ^ (i - e)- (7) 

Let us note, that the same arguments for the pseudoscalar wf imply its 
constant behavior at ^ — > ibl. 

Now we want to repeat these arguments for the analysis of the transverse 
distribution. To do so, let us define the mean values of the transverse quark 
distribution by the following matrix elements: 

where is the covariant derivative, acting on the one quark and trans- 
verse vector = (0, t, 0) is perpendicular to the hadron momentum q^ = 
iqo,0±,qz). The factor ^^(2n)l|" is introduced to (H) to take into account 
the integration over (p angle in the transverse plane: / (i0(cos0)^"/ J d(j) = 
(2n- l)!!/(2n)!!. 

We interpret the {k\) in this equation as a mean value of the quark 
perpendicular momentum. Of course it is different from the naive, gauge 
dependent definition like {0\d'yu'y5d'l_u\7i{q)) , because the physical transverse 
gluon is participant of this definition. However, the expression (^) is the 
only possible way to define the k\ in the gauge theory like QCD. We believe 
that such definition is the useful generalization of the transverse momentum 
conception for the interactive quark system. 
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The first few transverse moments defined above, can be calculated using 
some special QCD sum rules |]19|. To find the behavior (A;^") at large n we 



can repeat our previous duality arguments with the following result[|: 

Here we introduced the notation S'7r(-L) for the vr meson interval duality (do 
not confuse with corresponding notation for the longitudinal distribution 
5.(11)). 

The nice feature of (^ is its finiteness for arbitrary n. It means that the 
higher moments 



do exist. In this formula we introduced the nonperturbative wf tlj{k'j_,C,), 
normalized to one. Its moments are determined by the local matrix elements 
(H). The relations to Brodsky and Lepage notations bl{xi, k±) and to 
longitudinal distribution amplitude 0(^) introduced earlier, look as follow: 

^BL(xi,fc^) = ^^^^^(e,fc^), dkiij{ki,o = m, / d^m = xm 
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where f-^ = 133MeV. The existence of the arbitrary high moments (A^i") 
means that the nonperturbative wf, defined above, falls off at large trans- 
verse momentum k\ faster than any power function (exponential?). 

Another feature of the formula (^) is apparently the bad news, because 
it corresponds to a nonphysical 5(^1 — S',r(-L))-shape distribution. It is abso- 
lutely clear that the physical interactions spread out this behavior and will 
change this result .0 Indeed, as will be shown in the next section, the higher 
loop corrections give a desired factorial behavior on n. 

To anticipate the events we would like to formulate here the obtained 
result. Higher loop corrections will change the behavior (^ in the following 
way: 

s:^\±)^s:^\±){i + ...+c^{^rn\}, c„~i. (n) 

TT 

For the large enough n the perturbative corrections will play the dominant 
role and they actually define the large kj_ behavior of the nonperturbative wf, 
which turns out to be Gaussian. The relation (|^) converts this information 
into the fixing of the asymptotic behavior of wf at large kj_: 

(fcf ) = J dkldik^^ij{kl,C) ~ n! n^oo (12) 

=^ i^{kl ^ oo, x) ~ exp( ^ ), 

X 1 — x) 



^ Here and in what follows we ignore any mild (nonfactorial) n-dependence. As we 
discuss later, we will be able to find a factorial dependence, but not a prefactorial behavior. 
Thus, they will be consistently omitted. 

^ The situation here in some sense resembles the deep inelastic scattering. As is known 
the asymptotic form of the nucleon structure function at — > cxd is F2{x,q^ —> oo) ^ 
S{x). However, F2 differs strongly from the delta function and goes to this limit only at 
asymptotically high when the terms like 1 / In are neglected. 
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where we took into account the previous result [|T^ that ^^-^f^^ is the only 
allowed combination the wf can depend on. 

Few comments are in order. We identify the moments @ defined in 
QCD with the ones defined in the constituent quark model. We believe that 
such an identification is the useful working conjecture, which might be helpful 
in understanding the connection between QCD and the valence quark model. 

Now let me recall few results (from the constituent quark model) which 



support this conjecture. It is well known pO[ that the equal- time wave 
functions 

i'cui.^) ~ exp(-g^) (13) 

of the harmonic oscillator in the rest frame give a very reasonable description 
of static meson properties. Together with Brodsky-Huang-Lepage prescrip- 



tion [pll] connecting the equal -time and the light-cone wave functions by 



identification 



4a;(l — x)' '^^"'^'^ > '^'^^ 4x(l — x)^' 

one can reproduce the Gaussian behavior (|T2|) found exclusively from QCDQ. 
It means, first of all, that our identification of the moments (H) defined in 
QCD with the ones defined in quark model, is the reasonable conjecture Q. 

S.Perturbative series at large orders. 

This section is the technical one. Its purpose is the explanation of the 
formula (|ll]) which constraints the asymptotics of the wf in the transverse 
direction. 

The main idea, as before, is analysis of the appropriate correlation func- 
tions at large q^. We estimate the perturbative corrections at large n to find 
the dependence on external number n. Our main assumption corresponds 
to the very likely case that the tt meson gives a nonzero contribution to the 
dispersion integral for arbitrary n. Formally, it can be written in terms of 
dispersion relation (^. We assume that S',r in this relation is some finite, 
nonzero number. As soon as this assumption has been made, the standard 
calculation of the n dependence of correlation function convert this informa- 
tion into the knowledge about n dependence for the tt meson matrix elements 
(^^. This information is required for the fixing of the wf asymptotic be- 
havior. 



Before going on, let me pause here to recall Lipatov's method |22] of the 



calculation of large order behavior in a field theory (for recent review of this 

^More exactly, in this identification the combination with constituent effective quark 
mass k^+m? appeares. However, because we are interested in the asymptotic k\ behavior, 
we skipped these terms. 

^The same method can be apphed for the analysis of the asymptotical behavior of the 
nucleon wf which in obvious notations takes the form: 

^2 

1pnucleon{k\i ^ OO, Xi) - exp(- ^ -^) 
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technique, see [^). For great simplification we would like to demonstrate 
the main idea of the method in a scalar field theory. The generalization for 
the nonabelian gauge theories is straightforward. 

The main object of our analysis is the n-point correlation function, Z^'^\g), 
where n is the number of external fields (external currents) and g is the cou- 
pling constant in a scalar field theory. For g small, Z^{g) is given by an 
asymptotic expansion: 

Z^-\9) = J:zI''^9' (14) 

k 

We would like to calculate the large order behavior of this correlation function 
at /c — s> oo. The main Lipatov's idea is to present the coefficients zj^^ through 
a contour integral in the complex g— plane: 

where S{(j)) is the action of the scalar field theory and D(j) is the standard 
measure for the functional integral which defines the theory (We confined 
ourself in this formula the case n = only. It corresponds to analysis of the 
Grand Partition Function Z{g) itself. The dependence on n will be discussed 
a little bit later.). If the theory possess the classical instanton solution, 
then the calculation of the integral over g can be done through steepest 
descent method. This method is justified only for small g. But for k large, 
the integral over g is dominated by the small g contribution. In particular, 
for four dimensional (j)^ field theory the classical instanton solution has the 
property <^c/ ~ ^1 \f9-, [Hi and the classical action ~ l/s' is parametrically 
large. Thus, the semiclassical approximation is completely justified. 

In our case, when we consider the correlation function with external 
fields, and not the Grand Partition Function Zi^g^ itself, we have to substitute 
in the first approximation the classical solution 0^; in place of the external 
fields. More precisely, 

Zt^ -Jd^pJ -g^0(a;O0(x2)...0(x„)e-^W ~ (16) 

Jd<pJ -^e-^(*-VcK^i)0cz(x2)...0c/(a;„) ~ ^fcTTTf ~ + 

In this formula we took into account that the classical field depends on g as 
(pel ~ and the total number of external fields in the correlation function 

is equal n. It gives an additional factor g^ in denominator. The factorial 
dependence (|T6|) can be interpreted as the rapid growing of the number of 
Feynman graphs. In the spirit of the statistical approach to evaluation 
of the Feynman graphs one could interpret the same factorial dependence as 
the large multiplicity of order k of the virtual states. 

The semiclassical approximation we have used in the derivation (|16|) is 
justified as far as number (A; + |) 3> 1. Only in this case the integral over 
g is dominated by the small g contribution and instanton calculus can be 
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applied. It means that this method is still applicable when we calculate the 
dependence on the number of external fields n ^ 1. Thus, by using this 
technique we can find out the dependence on number n which is the main 
goal of the section. 

After this reverence to scalar field theory we come back to our analysis 



of the appropriate correlation functions in QCD. It can be shown |19| that, 
up to some unessential numerical factors, the analysis of transverse moments 
can be reduced to the consideration of the operators ~ q^o-^yGfj^^Yq. Thus, 
we arrive to the following formula for the n-dependence of the transverse 
correlation function: 

OA I ^'^-'"-^•"((G.J.)" (17) 



r\-/ 



Jo ' ~(fc + -)!,~n! 



where we took into account that for the classical instanton solution the field 
strength tensor is ^ 1/g. The classical action for instanton is Sd ~ 1/5'^- 
Besides that, we took into account that in order to provide the correct order 
in the perturbative expansion for n-th term ~ a", we need to take the k = ~ 
th order in the formula ([T7|). It gives the appropriate factor {g'^)^ = g"' from 
the interaction of the external gluons with themselves, see Fig.2. Finally, we 
arrive to the n! dependence in the formula (|T^, announced and intensively 
employed earlier (pi]). 

Few comments are in order. First of all, the similar analysis of the longitu- 
dinal correlator (^ does not lead to the n\ behavior found for the transverse 
correlation function. The technical reason for that is the selection of the 
special kinematical structure ~ (z^g^)", 2;^ = in (|^). An instanton, simply 
by kinematical reasons, does not contribute to this structure in the leading 
order. 

Now we wish to understand the difference, mentioned above, between lon- 
gitudinal and transverse correlators in terms of standard Feynman diagrams. 
As we discussed earlier, the factorial dependence in the field theory is the 
simple refiection of the rapid growing of the Feynman graphs. Thus, we have 
to be surprised not by the fact that such factorial dependence appeared, but 
rather, the fact that sometimes it does not show up. Therefore, let us try to 
understand why the extracting of the longitudinal structure suppresses the 
most of the graphs? 

To answer on this question, let us consider a typical graph with three- 
gluon vertex. Fig. 3 (exactly such kind of graphs, at least in the limit of 
large number of colors, gives the most important statistical contribution) 
and take the longitudinal projection for the external gluons ~ (A^z^)"^. We 
wish to extract the structure ~ (q^z^)'^ after completing the calculation. The 
direct computations demonstrate that the corresponding contribution to this 
structure is zero. The only nonzero result comes from the graphs like Fig. 4, 
where the external gluons interact with quarks, but not with themselves. It 
is clear, that the number of such diagrams does not demonstrate the factorial 
growth. 
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The explanation is simple: transverse and longitudinal structures are 
quite different because transverse gluons (in contrast to longitudinal ones) 
are physical degrees of freedom and can not be rotated away by a gauge 
transformation. Both methods presented above (instanton calculus and dia- 
grammatic approach) lead to the same conclusion and agree with each other. 

4. Nonlocal condensate and constituent mass of the quark. 

The nonlocal condensate (NLC) has been introduced into the theory 
(to my knowledge) within operator product expansion (OPE) technique by 
Shuryak more than 10 years ago |2^.| This object naturally arises in the 



analysis of the heavy-light quark system, qQ. Indeed, if we consider along 
with ||2^,[^ the correlation function {T{qQ{x), Qq{0)}), describing this sys- 



tem, we end up (in the limit Mq oo) with the object which is completely 
factorized from the heavy quark and which is called NLC: 



{T{qQ{x),Qq{0)}) ~ {q{x)P exp{ig / A^dx^)q{0)) + perturb, part. (18 



X 







All nontrivial, large distance physics of the system is hidden there. The NLC 
is gauge invariant, well defined object and thus, it deserves to be studied in 
a more detail. As we mentioned, this system may provide a definition of the 
constituent quark in QCD. We shall discuss this important point a little bit 
later, but now we would like to note that for the large Euclidean x the NLC 
is dominated by the lowest state contribution and thus [ p^ : 

{q{x)Pexp{ig f A^dx^)q{Q)) °° exp(-Av^), (19) 

•J 

where A is the lowest energy level of the mesons in the heavy quark effective 
theory (HQET) (see, e.g. review [p^ ]). 

Another interesting application of the NLC is the analysis of the hadron 
properties within QCD sum rules However, we shall not discuss 



any applications of the NLC in this letter. Instead we concentrate on the 
analysis of the NLC itself and its general properties. 

It is convenient to parametrize the coordinate dependence of the NLC by 
a Laplace-type representation pB| : 



{q{x)Pexp{tg A^dx^)q{0)) = {qq)S{x'), S{x') = e— f{s)ds. [20) 
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From its definition it is clear that the n— th moment of the spectral function 
f{s) can be expressed in terms of local operators with 2n derivatives: 

^i^rm I f{s)s-ds = ^M^.D.fy) (21) 

Thus, the entire information about NLC is hidden now in the properties of 
the spectral function f{s) or (which is the same) is coded by the behavior of 
the high dimensional condensates ~ (g(D^)^"g). 

^The term NLC itself appeared later, |26|. 
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What do we know about properties of the /(s)? First of all, we know its 
normalization and the first moment: 



/ f{s)ds = l, 
Jo 



, fis)sds = ^'Yf^ = ^'''y^r'^ = ^ ^ 0.2GeV^ (22) 
' ^ ^ 2{qq) Mm) 4 ' V ; 

where we introduced the standard notation (see e.g.0) rnl ~ O.SGeV^ for 
the mixed quark-gluon vacuum condensate of the dimension five. 

Besides that, we know the asymptotic of the NLC for large euclidean x 
(p!9|). The transformation (pO|) converts this information into the fixing of 
the asymptotic behavior of the spectral function at small s : /(s — > 0) ~ 
exp(-^). 

In addition to these, well established results, we have some more con- 
straints. First, we know the next, high dimensional condensate, discussed 
earlier [|19|)EB- Second, the large n- behavior for such condensates can be 



found by the method elaborated in the previous section. 

Let us start from the analysis of the next moment of the spectral func- 
tion f{s). The formula ( pTD gives the following relation between the second 
moment of f{s) and the vacuum condensate of dimension seven: 

-{qq) L f{s)s'ds={q{x,D,rq)= (23) 



4 '''' Jo 
—3 3 
x^{-^{qg'^(^^^uG^^(JXaGxaq) + —{qg^G^^Gf^^q)}, 

where we evaluated the condensate {qD^D^D^D^q) in a standard way, by 
reducing it to the operator which depends only on G^,^. [|Ty],|]2S|. To estimate 



the right hand side of the eq.(^), we introduce the coefficient of "nonfactor- 
izability" K for the condensates of dimension seven (one and the same for 
all operators). In particular, 

—K 

{qg^a^^G^^axaGxaq) = -^{m) {g'^G''^uG%) , 

{qg'G,M,,q) = j{qq){g'G%G;,). 

(In this formula i^' = 1 if the factorization prescription would work). The 
coefficient K was estimated numerically in [l^,[|29| and was found to be 
from 3 to 4. There is an additional, independent estimate which supports 
this value for K and which will be presented at the end of the section. 

Finally, we arrive to the following value for the second moment of the 
spectral function: 

POO 

f{s)s'ds = ^ 0.04 - 0.05GeV\ (24) 

where we use the standard value for the gluon condensate {^G'^^J) ~ 1.2 ■ 
lQ~'^GeV^, H and the estimated value for the coefficient K = 3 4. 
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Our next constraint comes from the analysis of the asymptotically distant 



terms in the formula ( pTI) . We use PCAC in order to relate the matrix 



elements and vacuum condensates: 

{^\dt^,{x^D^,f'-'u\iT) = Z^(0|g(x^D^)2"g|0), (25) 

Jtt 

By repeating the previous arguments we arrive to the following conclusion 
(see footnote 2 on page 4): 

{q{x,D,fy) ~ {x^r{q{D'rq) ~ {x^f {q{tga,,G .^Tq) ~ n\ (26) 

This result combined with the definition (^) gives the following constraint 
on the asymptotic behavior of the spectral function f{s): 

s''f{s)dsr^l. (27) 





Let us note the absence of the factorial dependence on the right hand side of 
(^). It means that f{s) falls off faster than any exponential function. We 
parametrize such a dependence at large s by step fuction: /(s) ~ 0(cr^ — s), 
where B(x)— function is defined in the standard way: B(x > 0) = 1, 0(x < 
0) = 0. 

Combining this result with the small s behavior in the simplest way, we 
arrive at the following Ansatz for the spectral density /(s): 

f{s) = N-e-^e{a^-s) (28) 

with three dimensional parameters: A^, A, a. 

The Fourier transformed function S{p) with Ansatz ([28| ) takes the follow- 
ing form: 



, ,~ ds v^-p? 1 -A^+P^ 



Few comments are in order. First of all, all formulae in this section, start- 
ing from (|T^), are derived in the Euclidean space (imaginary time), x^ < 
0, < 0. However we expect that the obtained analytical properties (pole, 
in particular) will not be changed after the continuation to Minkowski space 
is made. In this sense the formula ( pOf ) demonstrates that we have dealt with 
a "constituent quark" with mass A and width a. This determines the 
physical meaning of our phenomenological parameters introduced by Ansatz 

mn 

Our next comment is the qualitative remark that our Ansatz (p8| ) for 
the spectral density /(s) describes the well-localized function (at sq — o"^) 

^We expect, of course, that the constituent quark mass depends on momentum p^. 
We beheve that this dependence would come into the game if the corrections ^ 1/n in 



eq.(27)would be taken into account. In particular, more complex Ansatz for the spectral 
density /(s) f{^)J2k=o'^i^^^' 7^ does not change the asymptotical behavior, but 
leads to the momentum dependence of the constituent quark mass. Unfortunately, we can 
not reconstruct these, non-leading terms, by the described method. 
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almost for arbitrary parameters A, a. Thus, for any reasonable choice of 
these parameters we would expect that the following feature of the Ansatz 
m) is fulfilled: 



(/(rs/(s)c/s)2 (so; 



2 



1. (30) 



It should be considered as an independent check of the nonfactorizability of 
the mixed vacuum condensate of dimension seven with factor ^ 3 -^ 4. 

We are prepared for some numerical estimates now. The cr parameter is 
fixed by the mixed vacuum condensate (|22|) and is equal to cr^ ^ 0.3GeV^. 
The parameter A presumably should coincide with the constituent quark 
mass, i.e. A ~ 0.356*61^. We sketch the NLC on Fig. 5 where we use these 
numerical parameters. 

Our last remark is a very general observation that the origin for the 
Gaussian behavior of the light cone wf (|l^) (it corresponds to the harmonic 
oscillator wf in the quark model, see([T3|)) and the very specific behavior 
(^,^) of the vacuum nonlocal condensate is one and the same: It is the 
factorial dependence of the pion matrix elements ( p6D in the former case and 
the same factorial dependence of the vacuum condensates {q{D'^)"'q) ~ n\ in 
the latter one. It determines both these phenomena. 

5. Conclusion. 

Let me formulate the main results of this letter. 

• We formulated some constraints on the nonperturbative light cone wf. 
The analysis was based on the very general methods like dispersion relation 
and duality, and thus, the obtained results should be considered as direct 
consequence of QCD. We found the Gaussian asymptotics for the wf, which 
we believe, not accidentally coincides with the harmonic oscillator wf from 
valence quark model. 

• We reconstructed the NLC, the fundamental object of QCD, by analysing 
its few moments and the asymptotical behavior. 

We hope that this letter can motivate the future investigations in the 
following direction: How one can understand the valence quark model (which 
is the main tool in the hadronic physics for already three decades) in terms 
of QCD Lagrangian? . 
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Figure Captions. 

Fig.l The asymptotic behavior of the two -point correlation function at 
one loop level. 
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Fig. 2 The large order corrections to the correlation function which lead 
to n\ behavior. 

Fig. 3 A typical graph with three-gluon vertex which contributes to the 
transverse, but not longitudinal structure. 

Fig. 4 A typical graph which contributes to the longitudinal structure. 

Fig. 5 The NLC ^^(^^^'^^p^'^^/q^ Af,dx^)q(o)) _ ^^^-^ function of distance for 
parameter A ~ 0.35GeV. 
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